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CN ■ We present rest-frame i?-band galaxy luminosity function measurements for three 

O . different redshift ranges: 0.5 < z < 0.75, 0.75 < z < 1.0, and 1.0 < z < 1.5. Our 

measurements are based on photometric redshifts for ~ 3000 H-b&nd selected galaxies 



in 



with apparent magnitudes 17 < H < 20 from the Las Campanas Infrared Survey. We 
show that our photometric redshifts are accurate with an RMS dispersion between the 
photometric and spectroscopic redshifts of cr z /(l + z) ~ 0.08. Using galaxies identified 
f^. ' in the Hubble Deep Field South and Chandra Deep Field South regions, we find, re- 

spectively, that (7.3 ± 0.2) % and (16.7 ± 0.4) % of the H < 20 galaxies are at z > 1. 
■ We first demonstrate that the systematic uncertainty inherent in the luminosity func- 

tion measurements due to uncertainties in photometric redshifts is non-negligible and 
therefore must be accounted for. We then develop a technique to correct for this sys- 
tematic error by incorporating the redshift error functions of individual galaxies in the 
Q-f luminosity function analysis. The redshift error functions account for the non-gaussian 

characteristics of photometric redshift uncertainties. They are the products of a con- 
volution between the corresponding redshift likelihood functions of individual galaxies 
and a Gaussian distribution function that characterizes template-mismatch variance. 
We demonstrate, based on a Monte Carlo simulation, that we are able to completely 
recover the bright end of the intrinsic galaxy luminosity function using this technique. 
Finally, we calculate the luminosity function separately for the total ii-band selected 
sample and for a sub-sample of early-type galaxies that have a best-fit spectral type of 
E/S0 or Sab from the photometric redshift analysis. The primary results of this analysis 
are: (1) the galaxy luminosity functions are consistent with a Schechter (1976) form; 
(2) the evolution of the co-moving luminosity density £ of ff-band selected galaxies is 
characterized by Alog^ /Alog(l + z) = 0.6 ± 0.1 at rest-frame 6800 A; and (3) £r for 
color-selected early-type galaxies exhibits moderate evolution from z ~ 1.5 to z ~ 0.3. 
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Specifically, £r for these red galaxies brighter than 1.0 (1.6) could decrease with in- 
creasing redshift by at most a factor of three (six) over this redshift range after removing 
possible stellar fading according to the most extreme stellar evolution scenario. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: luminosity 
function — surveys 

1. Introduction 

The galaxy luminosity function represents the product of a series of physical processes that 
govern galaxy evolution, including star formation and the recycling of metals through various 
heating and cooling processes. An accurate estimate of the galaxy luminosity function as a function 
of time therefore bears significantly on our understanding of how structures in the universe form 
and evolve. Independent measurements of the optical galaxy luminosity function based on various 
deep redshift surveys between z ~ and z ~ 0.75 have yielded consistent results (Lilly et al. 1995; 
Ellis et al. 1996; Marzke & da Costa 1997; Bromley et al. 1998; Marzke et al. 1998; Lin et al. 
1999; Madgwick et al. 2001). For example, the galaxy luminosity function is well represented by 
a Schechter function (Schechter 1976) with a steeper faint-end slope for galaxies of bluer optical 
colors at all epochs that have been studied. Comparisons of these measurements at different 
redshifts further show that these blue galaxies exhibit a stronger evolution than red galaxies both 
in their space density and in their rest-frame optical luminosity. 

It is difficult, however, to apply these results to fully constrain galaxy formation models. While 
galaxy luminosities measured in the optical are sensitive to the transient star formation rate and 
therefore are a good tracer of the instantaneous star formation activity, they are not represen- 
tative of the total stellar mass that has been assembled in the galaxies over time. Additional 
caveats exist in all studies that attempt to draw conclusions based on comparisons of luminos- 
ity function measurements obtained using different galaxy samples. For example, comparison of 
morphologically-selected and color-selected galaxy samples is inadequate, because even morpholog- 
ically distinct galaxies can exhibit similar colors. Furthermore, it is impossible to unambiguously 
identify the same galaxy population at different epochs using the observed colors because the un- 
derlying stellar evolution sequence is unknown. Finally, other selection biases such as the cosmic 
surface brightness dimming effect (preventing accurate morphological classifications at high red- 
shift), bandpass selection effect (galaxies selected in the same observed- frame bandpass may have 
different rest-frame properties), and the presence of dust (in which a young, star- forming galaxy 
masquerades as an old, quiescent one) further complicate the interpretations of comparisons of 
different galaxy surveys. 

We have initiated the Las Campanas Infrared (LCIR) survey (Marzke et al. 1999; McCarthy 
et al. 2001; Chen et al. 2002; Firth et al. 2002) to study the history of mass assembly based on 
a uniform sample of near-infrared selected galaxies. The rest-frame near-infrared light is a good 
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proxy for the total stellar mass, independent of age or galaxy type (Chariot 1996). Measurements 
of the galaxy luminosity function at rest-frame near-infrared wavelengths at different epochs may, 
therefore, be adopted to estimate the stellar mass density evolution. A robust measurement of the 
rest-frame near-infrared luminosity function at z ~ 1 requires a galaxy survey conducted in the 
K band. We have so far completed the LCIR H-b&nd survey and a complete .FT-band survey is 
underway. Here we present measurements of the rest-frame .R-band luminosity functions in three 
redshift intervals over the range z = 0.5-1.5 for the H-hand selected galaxies using photometric 
redshifts. At z ~ 1.2, the observed-frame if-band approximately corresponds to rest-frame i?-band, 
and therefore the calculation does not depend sensitively on the templates adopted to determine the 
^-correction. The primary objectives of the analysis are (1) to obtain measurements of the galaxy 
luminosity function at z > 1, where different galaxy formation scenarios begin to show distinct 
predictions in the space densities and masses of evolved galaxies, and (2) to study galaxy evolution 
at redshifts between z ~ 0.5 and z ~ 1.5 based on a uniform sample of near- infrared selected 
galaxies, thus extending the redshift range of existing deep redshift surveys to beyond z = 1. 

We first conducted a photometric redshift survey in the LCIRS galaxy sample to identify a 
complete sample of galaxies at z > 0.5. The photometric redshift survey in two of the survey fields 
identified > 3000 galaxies at z > 0.5, yielding a statistically significant sample of intermediate- 
redshift galaxies for understanding the galaxy population at redshifts beyond z = 1. The sample 
size is comparable to those of wide-field spectroscopic surveys carried out at z < 0.5 such as 
the Canadian Network for Observational Cosmology Field Galaxy Redshift Survey (CNOC2; e.g. 
Lin et al. 1999). The photometric redshift techniques determine galaxy redshifts based on the 
presence/absence of broad-band spectral discontinuities, rather than narrow-band emission or ab- 
sorption line features. They are therefore of particular value for the identification of galaxies at 
z > 0.75. Complete spectroscopic surveys of faint galaxies in this redshift range are extremely 
challenging, as prominent spectral line features in galaxy spectra are redshifted to near-infrared 
wavelengths and are difficult to detect due to bright, contaminating sky lines and detector limita- 
tions. It has been demonstrated over the past several years that distant galaxies may be robustly 
identified using photometric redshift techniques that incorporate optical and near-infrared broad- 
band photometry (Connolly et al. 1997; Lanzetta, Fernandez-Soto, & Yahil 1998; Fernandez-Soto, 
Lanzetta, & Yahil 1999; Martini 2001; Fernandez-Soto et al. 2001; Rudnick et al. 2001). In this 
paper, we show that the photometric redshift measurements of the galaxies identified in the LCIR 
survey are accurate, with an RMS dispersion between the photometric and spectroscopic redshifts 
of a z /(l + z)& 0.08. 

Next, we developed a maximum likelihood method to measure galaxy luminosity function 
using photometric redshifts, which explicitly accounts for photometric redshift uncertainties of 
individual galaxies. Photometric redshift uncertainties together with a steep slope of the galaxy 
luminosity function would flatten the observed luminosity function and yield systematic biases. 
Uncertainties in photometric redshifts must therefore be accounted for in the luminosity function 
analysis. We demonstrate that systematic uncertainties in the derived galaxy luminosity function 
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are significantly reduced when redshift error functions are taken into account. We show that (1) the 
galaxy luminosity function is consistent with a Schechter (1976) form; (2) the evolution of the co- 
moving luminosity density of the H-h&nd selected galaxy is characterized by Alog£ /Alog(l + z) = 
0.6 ± 0.1 at rest-frame 6800 A; (3) galaxies with spectral energy distributions best characterized 
as E/SO or Sab appear to have faded by ps 0.5 mag from z ~ 1 to z ~ 0.5, consistent with the 
expected luminosity evolution from an exponentially declining star formation rate model; and (4) 
only moderate evolution in the space density and in the rest-frame, co-moving i?-band luminosity 
density is found for the bright, color-selected early-type galaxies. 

We adopt the following cosmology: VLyi = 0.3 and Vt\ = 0.7 with a dimensionless Hubble 
constant h = i?o/(100 km s _1 Mpc -1 ) throughout the paper. 

2. The Las Campanas Infrared Survey 

The LCIR survey is a deep, wide-field near- infrared and optical imaging survey designed to 
identify a large number of galaxies dominated by an old stellar population at z > 1, while securing 
a uniform sample of galaxies of all types to z ~ 2 based on their broad-band optical and near- 
infrared colors. We have completed the H-band survey, covering one square degree of sky to a 
mean 5<r detection limit in a four arcsecond diameter aperture of H ~ 20.8. Initial catalogs of > 
54,000 galaxies identified over 1400 arcmin 2 together with results from a comprehensive study of 
the survey incompleteness in each field are presented in Chen et al. (2002). The analysis presented 
in this paper is based on galaxies identified in two of the LCIR fields: the Hubble Deep Field South 
(HDFS) and Chandra Deep Field South (CDFS). To summarize, the galaxy sample contains ~ 
6,700 ii-band selected galaxies over 847 arcmin 2 in the HDFS region with complementary optical 
U, B, V, R, and I colors, and ~ 7, 400 ii-band selected galaxies over 561 arcmin 2 in the CDFS 
region with complementary optical V, R, I, and z' colors. 

To study the statistical properties of galaxies at z > 1 using the LCIR survey sample, we need to 
first determine redshifts for the galaxies. The photometric redshift techniques that determine galaxy 
redshifts based on their broad-band spectral energy distribution (SED) offer an efficient means to 
obtain robust redshift measurements together with well-understood measurement uncertainties for 
all galaxies in the survey. 

2.1. Photometric Redshift Likelihood Analysis 

We adopted the photometric redshift technique originally developed by Lanzetta and collabo- 
rators (e.g. Lanzetta, Yahil, & Fernandez-Soto 1996), and modified the program for its application 
to ground-based imaging data. The technique determines galaxy redshifts by comparing the SED 
of a galaxy, which is established based on a compilation of broad-band photometric measurements, 
with a grid of spectrophotometric templates evaluated at different redshifts. It employs a maximum 
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likelihood method, in which the likelihood of a galaxy matching a given spectral type T at redshift 
z based on n photometric measurements is written as 



C z (z,T) = J|exp 



i=l 



fi-AFi(z,T) 



1 2 



(1) 



where fi and Fi are respectively the measured and model galaxy fluxes in bandpass i, A is the 
normalization of the template, and Oi is the measurement uncertainty of fi. In cases where the 
galaxy is not observed in m of the n bandpasses to the detection limits of the images, equation (1) 
is modified as 



C z (z,T) = exp 



fi-AFi(z,T) 



1 2 



Oi 



l J max 

n / 

U = l •'/mm 



f'-AFi(z,T) 



Oi 



(2) 

where / m j n and / max are the minimum and maximum possible observed fluxes of the galaxy. We 
have chosen / m j n = and / max = 2<7j. The likelihood function C z (z,T) is maximized with respect 
to the type T at different redshifts to form C z (z), which is then maximized with respect to the 
redshift z to determine a best-fit photometric redshift z p hot for the galaxy. 



2.2. Spectral Energy Distribution Templates 

We included in our photometric redshift likelihood analysis six galaxy templates, a suite of 
stellar templates, and a QSO template in order to effectively identify the large number of stars and 
potential QSO/AGN candidates that appear in wide- field surveys. Every template covers a spectral 
range from A = 300 to 25,000 A. The galaxy templates were adopted from Fernandez-Soto et al. 
(1999) and Yahata et al. (2000) and spanned a range of spectral types — from elliptical or SO (E/S0), 
Sab, Scd, irregular (Irr), to starburst (SB). A total of 150 stellar templates were compiled from the 
literature, spanning a wide range of stellar types — from early-type O and B stars to late-type L and 
T dwarfs (Oppenheimer et al. 1998; Pickles 1998; Fan et al. 2000; Leggett et al. 2000). The QSO 
template was formed by first adopting the composite QSO spectrum presented in Francis et al. 
(1991) and extrapolating the spectrum to infrared wavelengths using a simple power law. Because 
the Francis spectrum was derived based in part on a large number of QSOs at z > 2 and was not 
corrected for intervening absorption at wavelengths shortward of the redshifted Lya emission line, 
we replace the continua at A < 1450 A with the composite spectrum presented in Zheng et al. 
(1997) based on ~ 100 spectra of QSOs at z < 1.6 obtained with the Faint Object Spectrograph 
aboard the Hubble Space Telescope. Our approach to separate stars from galaxies is described in 
Chen et al. (2002). The number- magnitude diagram of the stars identified in the HDFS and CDFS 
fields, in comparison to that of the galaxies, is presented in Chen et al. as well. 
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2.3. The Photometric Redshift Survey 

We performed the redshift likelihood analysis and determined photometric redshifts for all 
galaxies identified in the LCIRsurvey using the available near-infrared and optical broad-band pho- 
tometry. As discussed in § 2.4, uncertainties in galaxy photometry contribute to the uncertainties in 
photometric redshifts. Because galaxies vary in size, fixed-aperture photometry would include ex- 
cess noise from sky for smaller objects and exclude flux from the outskirts of bigger objects, both of 
which would yield a reduced signal-to-noise (S/N) ratio in the photometric measurements. Fixed- 
aperture photometry is, therefore, not optimal for establishing the observed SEDs of individual 
galaxies to be compared with model templates in the photometric redshift analysis. 

It is clear that to improve the accuracy of photometric redshifts, we need to first achieve an 
optimal S/N ratio for the photometric measurements of individual galaxies in all bandpasses, where 
by optimal measurements we mean those that do not include excess noise or exclude object fluxes. 
To reach the goal, we obtained a new set of photometric measurements using a varying-size aperture 
for different objects, which is different from the fix-aperture photometry published in Chen et al. 
(2000). The individual object apertures were adopted from the object segmentation determined by 
applying SExtractor (Bertin & Arnout 1996) to the "white light" image of each field. As described 
in Chen et al. (2002), the white light image is an unweighted sum of all the registered optical images 
that had been scaled to unit exposure time and filter throughput. Because of the improved signal, 
particularly in the low surface brightness regions, this image is more sensitive than any individual 
optical images for SExtractor to identify the true extent of each object at a given surface brightness 
detection threshold. We measured, for each galaxy, the total fluxes and associated uncertainties in 
different bandpasses using the same aperture determined in the white light image by SExtractor. 
This procedure allowed us to recover the total fluxes of each object in different bandpasses, without 
significantly compromising the precision of the photometric measurements. 

Note that the main objective of the analysis here is to obtain accurate and precise broad- 
band colors. Accurate measurements of the total fluxes of individual galaxies are necessary for a 
robust estimate of the galaxy luminosity function, but not for the purpose of obtaining an accu- 
rate photometric redshift. We adopted the optimal photometric measurements for the photometric 
redshift analysis and the 4" -aperture photometric measurements for the luminosity function anal- 
ysis. We present in Figure 1 the median residuals and the corresponding 1 a scatters between 
the "optimal" and 4" -aperture photometric measurements versus magnitude for the i/-band se- 
lected objects in the CDFS region in the V, R, I, and H bands. It was shown in Chen et al. 
(2002) that the 4" -aperture photometric measurements closely resembled the total magnitudes of 
individual galaxies — only < 0.1 mag aperture correction was necessary. Figure 1 shows that the 
optimal photometry is consistent with the fixed aperture photometry in all bandpasses to within 
measurement uncertainties, although it systematically underestimates the total fluxes at faint lim- 
its. We note that the systematic biases of the optimal photometry at faint magnitudes do not affect 
our photometric redshift estimates, because they are almost entirely independent of the adopted 
bandpasses. 
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The products of the photometric redshift analysis for each galaxy include: the photometric 
redshift measurement, the redshift likelihood function, and the best-fit spectral template. Figure 
2 shows examples of the results in pairs of panels for four HDFS galaxies that have been spec- 
troscopically identified at different redshifts. The top panel of each pair shows the photometric 
measurements and uncertainties as solid points with error bars. The solid curves represent the 
best-fit spectral templates with the spectral type indicated in the upper-left corner. The open 
squares indicate the predicted fluxes in individual bandpasses, estimated by integrating the best-fit 
template over the corresponding system transmission functions. The dashed lines indicate the zero 
flux level. The bottom panel of each pair shows the redshift likelihood functions. The estimated 
photometric redshift (z p hot) and spectroscopic redshift (z S pcc) are indicated in the lower-right corner 
of each of the top panels. In all four cases, we have successfully determined the galaxy redshifts to 
within the estimated uncertainties. 

We studied the redshift distribution of the .ff-band selected sample. The left panel of Fig- 
ure 3 shows the redshift distribution of the if-band selected HDFS galaxies with H < 20. The 
photometric redshifts were obtained based on the available U, B, V, R, I, and H photometric 
measurements. The hatched region indicates the redshift distribution of galaxies with I — H > 3. 
The I — H > 3 color criterion is predicted by various evolutionary scenarios to identify evolved 
galaxies at z > 1 (Chen et al. 2002). The curves indicate the predicted redshift distributions of 
the total i^-selected sample and the red population, which were calculated using the model that 
best fit the number-magnitude relation presented in Chen et al. (2002). Specifically, we calculated 
the predicted redshift distribution of the total ff-band selected galaxies by adopting the i"C-band 
luminosity function obtained by Gardner et al. (1997), assuming a rest-frame color of H — K ~ 0.2 
(e.g. Mobasher et al. 1986) and a no-evolution scenario. The predicted redshift distribution of the 
red galaxies was estimated based on the observed population fraction of elliptical galaxies in the 
local universe (e.g. Binggeli, Sandage, & Tammann 1988). We first adopted the same luminosity 
function and scaled it according to M*(ied) = — 1.2 and 0*(red) = 0.2 0* for identifying the red 
population. Next, we assumed an exponentially declining star formation rate history with a 1-Gyr 
e- folding time for galaxies formed at Zf = 30 to model the pure luminosity evolution and added 
random photometric noise of ~ 0.25 mag in the model to account for the uncertainty of the I — H 
color selection criterion in the data. Finally, we calculated the integrated galaxy number counts 
versus redshift. Figure 3 shows that these models agree well with the observation. 

In the right panel of Figure 3, we present the cumulative redshift distributions of the H < 20 
galaxies in the HDFS (open) and CDFS (hatched histogram) regions. We exclude the CDFS 
galaxies identified at z < 0.75 from our analysis as they lack U- and i?-band data, which causes 
photometric redshifts to be less accurate in this range. We find that (7.3 ± 0.2) % (HDFS) and 
(16.7 ± 0.4) % (CDFS) of the ff-band selected galaxies with H < 20 are at z > 1. The errors were 
estimated using a Monte Carlo simulation, which takes into account photometric uncertainties 
and photometric redshift uncertainties, the subject of the next section. The difference between 
the two fields in the number density of z > 1 galaxies is likely due to the fluctuations in large- 
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scale structures. In particular, our survey is sensitive to galaxies at the bright end of the galaxy 
luminosity function, which are believed to be more strongly clustered. 

2.4. Uncertainties in Photometric Redshift Measurements 

In order for any statistical analysis of the intrinsic properties of galaxies identified using photo- 
metric redshift techniques to be plausible, we must first quantify the uncertainties. As discussed in 
detail in Lanzetta et al. (1998) and Fernandez-Soto et al. (2001, 2002), uncertainties in photometric 
redshift measurements are due to two sources: (1) for bright galaxies, the dominant uncertainty is 
due to the finite number of templates employed in the analysis; this is known as template-mismatch 
variance, and (2) for faint galaxies, an additional uncertainty is due to photometric measurement 
uncertainties. Fernandez-Soto et al. (2001) demonstrated that uncertainties in photometric red- 
shifts due to template mismatch is well characterized by a Gaussian distribution function with the 
1-a width equal to the RMS residuals between photometric and spectroscopic redshifts. In addition, 
equation (1) explicitly shows that in the absence of template- mismatch variance the uncertainty 
of the redshift likelihood function depends entirely upon the accuracy and precision of the flux 
measurements. Namely, the smaller the cr^'s are, the more precise the best-fit photometric redshift 
measurement. Uncertainties in photometric redshifts due to photometric errors are, therefore, well 
represented by the redshift likelihood functions of individual objects. 

Photometric redshift uncertainties can be quantified by explicitly taking into account the red- 
shift likelihood functions without adopting a parametric form. The error function of the photometric 
redshift measurement of galaxy i is 

p i {z-z i ;z i )= / j^C^z' -z l -z i )-^ V [-{z-z'f/2a 2 z \dz'. (3) 

Jo a z \J 2ir 

Equation (3) is a convolution of the redshift likelihood function C\{z) and a Gaussian kernel of 
width a z = a z (l + z-i) that characterizes the template-mismatch variance, where a z is the RMS 
residuals between spectroscopic and photometric redshifts at zero redshift. 

To assess the accuracy and reliability of our photometric redshifts, we first compared our 
photometric redshifts with known spectroscopic redshifts. There are over 100 objects in the HDFS 
region (Glazebrook 1998; Tresse et al. 1999; Palunas et al. 2000; Cristiani et al. 2000; Dennefeld 
2002). Excluding the ambiguous spectroscopic redshift identifications indicated by these authors, 
we found 67 objects that have both reliable spectroscopic and photometric redshift measurements 
at z < 2.5, 13 of which are stars. There are 162 known spectroscopic redshifts in the CDFS region 
obtained by ourselves (McCarthy et al. 2002, in preparation). Excluding ambiguous photometric 
redshifts (objects with uncertain redshift likelihood distributions), we found 140 objects with both 
reliable spectroscopic and photometric redshifts. 

Figure 4 shows the comparisons of photometric (z p hot) and spectroscopic (z spec ) redshifts for 
the HDFS and CDFS objects in the left panel. Solid points represent objects identified in the HDFS 
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region and open circles represent those in the CDFS region. The solid line indicates z p hot = z spec . 
The distributions of redshift residuals are shown in the right panel of Figure 4 for the HDFS (open 
histogram) and CDFS (shaded histogram) objects. We find an RMS dispersion of a z /(l + z) ~ 0.08 
for the 67 HDFS objects at z < 1. Only one of the 13 stars was misidentified as a galaxy at z p hot = 
0.09. Note that the mean redshift residuals are nearly zero at all redshifts that have been tested 
here, supporting that there is no systematic bias in our photometric redshifts. It demonstrates that 
the photometric redshifts obtained using UBVRIH photometry are accurate at all redshifts below 
z ~ 1. The comparison based on the 34 CDFS objects at z > 0.75 shows <Xz/(l + z) ~ 0.08. It 
demonstrates that the photometric redshifts obtained using VRIz'H photometry are accurate at 
redshifts beyond z ~ 0.75, although they become unreliable at lower redshifts due to the lack of U 
and B photometry. 

There are only a small number of known spectroscopic redshifts at 1 < z < 2.5 in both the 
HDFS and CDFS regions, but it is unlikely that the accuracy of the photometric redshifts at 
1 < z < 2.5 would be much worse than those at lower redshifts. It has been demonstrated by 
different groups that photometric redshifts are accurate to within a z /(l + z) < 0.1 at 1 < z < 1.5 
for all galaxies that have been tested so far (Connolly et al. 1997; Fernandez-Soto et al. 2001; 
Rudnick et al. 2001). In particular, Fernandez-Soto et al. (2001) showed that using the same set 
of SED templates their photometric redshifts are consistent with spectroscopic redshifts for all 19 
galaxies that have known spectroscopic redshifts at 1 < z < 1.5 in the HDF north. This supports 
that the number of templates in our photometric redshift analysis is sufficient to minimize template 
mismatch uncertainties. A more critical examination of the accuracy of our photometric redshifts 
in this redshift range relies on the availability of more spectroscopic measurements. 

Next, we performed a series of simulations designed to evaluate the uncertainties in the pho- 
tometric redshifts due to purely photometric uncertainties. We first generated a set of photometric 
measurements in each bandpass for a given If -band magnitude and input redshift z- ul based on an 
E/S0 template. We then perturbed the photometric measurements in individual bandpasses within 
the 1(7 photometric uncertainties measured for the individual optical and H images. Next, we 
determined a best-fit redshift z out using the photometric redshift technique. Finally, we repeated 
this procedure 1000 times and calculated the histogram of the residuals between z- m and z ou t- We 
repeated the entire process for the Sab, Scd, and Irr spectral templates to assess the performance 
of the photometric redshift technique for galaxies of different spectral shape. Figure 5 shows an 
example of the simulation results based on the filter combination and image sensitivities of the 
HDFS data. The left panels show the histograms of the residuals between the input galaxy redshift 
Zi n and the best-fit photometric redshift z out and the right panels show the cumulative fractions 
versus absolute redshift residuals, = |^out — ^in 

I . The results indicate that while we could accu- 
rately identify an Sab-Irr galaxy of H = 20 at z = 1.2, there is an ~ 10% chance that we would 
misidentify an E/S0 galaxy of H = 20 at z = 1.2 as a starburst galaxy at z ~ 5. The bimodal 
distribution of the residuals is due to the confusion between a 4000- A spectral discontinuity at low 
redshifts and a Lya continuum break at high redshifts. 
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We repeated the simulations for a grid of redshift and .ff-band magnitude based on the filter 
combination and image sensitivities of both HDFS and CDFS data. The results are presented in 
Figure 6, each panel of which now shows the sum of the residual histograms over all four spectral 
types. Figure 6 indicates that, even in the absence of template-mismatch variance, photometric 
redshifts become progressively more uncertain at fainter apparent magnitudes and higher redshifts, 
as indicated by the broader distribution of the residuals. The bimodal distribution that appears, for 
example, for an E/SO galaxy of H = 20 at z = 1.2 in the HDFS region does not exist for a galaxy of 
the same brightness at the same redshift range in the CDFS region, because the CDFS images reach 
to fainter magnitude limits. It is also clear, however, that there exists a non-negligible probability 
that galaxies at z < 0.4 in the CDFS region would be misidentified as galaxies at z ~ 3 even at 
bright magnitudes. The misidentifications occur primarily in cases when the input spectral template 
was an early-type galaxy. Because of the lack of U- and S-band photometric measurements for 
galaxies in the CDFS region, the observed SEDs are poorly constrained at rest-frame wavelengths 
shortward of the 4000-A break. The ambiguity disappears for the CDFS galaxies at z > 0.8, where 
the available F-band photometry begins to impose a strong constraint on the shape of the observed 
SEDs. 

In Figure 7, we present the cumulative fractions of the simulated photometric redshifts versus 
absolute redshift residuals for the four z = 1.2 cases presented in Figure 6. The left panels show the 
simulation results for the HDFS galaxies and the right panels for the CDFS galaxies. At z = 1.2, 
it is clear that photometric redshifts are accurate to within Az = 0.1 for the CDFS galaxies of 
H < 21, and that they become very uncertain for the HDFS galaxies at H ~ 21 because of 
their relatively larger photometric uncertainties. Specifically, there is a ~ 30% chance that the 
photometric redshifts of these galaxies are off by more than Az = 0.1. 

In summary, we conclude that the photometric redshift measurements of the HDFS objects 
based on U, B, V, R, I, and H photometry are robust for galaxies with H < 19 at all redshifts 
z < 1. The photometric redshift measurements of the CDFS objects based on V, R, I, z', and 
H photometry are robust for galaxies of H < 20 at all redshifts 0.75 < z < 1. We exclude 
CDFS objects with photometric redshifts z < 0.75 from subsequent analysis because of the lack 
of photometric measurements in the observed U and B bands (see also Connolly et al. 1995). 
The robustness of photometric redshifts at z > 1 remains to be examined more critically, because 
only a small number of known spectroscopic redshifts are available in our survey fields. However, 
consistent redshift estimates of the HDF galaxies at 1 < z < 1.5 support the photometric redshift 
technique using the six templates in our analysis for this redshift interval (e.g. Fernandez-Soto et al. 
2001). We note that redshift uncertainties due to dust reddening is not explicitly accounted for in 
our analysis. We consider the effect of dust as part of sources that lead to template mismatch. To 
further improve the reliability of the photometric redshift measurements and to reduce the number 
of low-redshift galaxies that are mis-identified as star-forming galaxies at high redshifts due to 
photometric uncertainties, we restricted the redshift range of the photometric redshift likelihood 
analysis to z < 2.5. This is appropriate because the fraction of H < 20 galaxies at z > 2.5 is 
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~ 0.1% (Lanzetta et al. 1999; Cohen et al. 2000; Shapley et al. 2001). 

3. Rest-frame i?-band Galaxy Luminosity Function 

Here we present the rest-frame .R-band luminosity function for three redshift ranges: 0.5 < z < 
0.75, 0.75 < z < 1.0, and 1.0 < z < 1.5. We emphasize that the results are based on photometric 
redshifts for ~ 3000 -ff-band selected galaxies with apparent magnitude 17 < H < 20. We calculated 
the luminosity function separately for the total .ff-band selected sample and for a sub-sample of 
early- type galaxies that have a best-fit spectral type of E/S0 or Sab in the photometric redshift 
analysis. The primary objectives are (1) to obtain an estimate of the galaxy luminosity function at 
z ~ 1 in order to study the statistical properties of galaxies at this epoch, and (2) to study galaxy 
evolution at redshifts 0.5 < z < 1.5 based on a uniform sample of near-infrared selected galaxies; 
this extends the redshift range of existing deep redshift surveys to beyond z = 0.75. 

The absolute i?-band magnitude Mr of a galaxy with an apparent magnitude H at redshift z 

is 

M R (z) = H- 25.0 - 5 x log - k'(z) + 2.5 x log(l + z), (4) 

where Dl{z) is the luminosity distance to the galaxy in units of Mpc and k'(z) is the /c-correction 
term to account for the color difference between the observed-frame and corresponding rest-frame 
bandpasses. The /c-correction term for each galaxy was calculated using the best-fit spectral tem- 
plate from the photometric redshift analysis. For galaxies at z ~ 1.2, the observed- frame H band 
approximately corresponds to rest-frame R, so the calculation does not depend sensitively upon the 
adopted templates to determine the /c-correction. For galaxies at z < 1, the estimated /c-correction 
ranges from between -0.34 (SB) and -0.9 (E/S0) at z = 0.5 to between -0.12 (SB) and -0.43 
(E/S0) at z = 1.0. 

Several methods to calculate the luminosity function have been proposed over the years. We 
examine three of these: (1) the 1/V max approach (Schmidt 1968; Felten 1976), (2) the Sandage- 
Tammann-Yahil (STY) approach (Sandage, Tammann, & Yahil 1979), and (3) the stepwise maxi- 
mum likelihood (SWML) approach (Efstathiou, Ellis, & Peterson 1988) in the context of a photometric- 
redshift determination of the luminosity function. A foreseeable systematic error in the derived 
luminosity function using photometric redshifts is to overestimate the number density of galax- 
ies at the bright end, producing a flattened, deduced galaxy luminosity function relative to the 
intrinsic luminosity function. This is due to the relatively large uncertainties associated with pho- 
tometric redshifts and the steep slope at the bright end of the galaxy luminosity function. There 
are more intrinsically fainter galaxies scattered into the brighter end of the luminosity function 
than intrinsically brighter galaxies scattered into the fainter end, thus yielding an apparently flat- 
ter luminosity function. To obtain robust estimates of the "intrinsic" galaxy luminosity function 
using photometric redshifts, it is therefore important to account for these systematic biases. In the 
following sections we present a maximum likelihood analysis that explicitly takes into account the 
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empirical redshift error functions of individual objects in the luminosity function calculation, and 
demonstrate, via a Monte Carlo simulation, that our approach allows us to reduce the systematic 
errors and to completely recover the bright end of the galaxy luminosity function. We summarize 
these three different luminosity function estimators, describe how we quantify and incorporate the 
photometric redshift uncertainties. 



3.1. The 1/Vmax Approach 

The 1/Vmax approach was first proposed by Schmidt (1968) and modified by Felten (1976). 
The galaxy luminosity function </>(M) is related to the number of galaxies per co-moving volume 
expected in a survey defined by the limiting apparent magnitude range m m \ n < m < m max and the 
redshift range z m \ n < z < z max according to 

An = 0(M)(M max - M min ) = J2^r, for M min < Mi < M max and z mhl < z { < z max . (5) 

i 

The absolute magnitude range [M mm , M max ] corresponds to the apparent magnitude range imposed 
by the survey at the galaxy redshift z^. The co- moving volume Vi is the maximum accessible volume 
of the survey for a galaxy of Mj, 

v = [*°» c0i D\ dz 
1 H ° (! + z ) 2 V^m(1 + zf - (fi M + n A - 1)(1 + zf + Q A ' 

where c is the speed of light, 0j is the angular area covered by the survey that is sensitive enough 
to detect the galaxy Mj at Zi, and z' min and z' max are the minimum and maximum redshifts that 
satisfy both the limiting magnitude and redshift ranges of the survey. The angular area A{ for each 
galaxy identified in the LCIR survey is determined from the results of the incompleteness analysis 
presented in Chen et al. (2002). 

The 1/Vmax approach is essentially a maximum likelihood method that does not assume any 
parametric form for estimating the galaxy luminosity function (see also § 3.3 below). It has been 
shown, however, that large-scale fluctuations in the galaxy number density introduce systematic 
biases in the galaxy luminosity function obtained with this approach for a magnitude-limited sample 
(de Lapparent, Geller, & Huchra 1989). We show in § 3.4 that it is difficult to recover the intrinsic 
galaxy luminosity function with this approach for a galaxy sample identified in a photometric 
redshift survey. 



3.2. The Sandage-Tammann-Yahil Approach 



Previous studies have demonstrated that the galaxy luminosity function may be well repre- 
sented by a Schechter (1976) function, 

$(L; L*) = 0* • {LjUY ■ exp(-L/L*) (7) 
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or 

$(M;M*) = (0.41nl0) • fa ■ W 0A W.-M)(i+ a ) . exp (_ 10 o.4 (m.-m)^ (8) 

To accurately determine the faint-end slope a and M* for the Schechter luminosity function in the 
presence of large-scale density fluctuations, Sandage, Tammann, & Yahil (1979) first introduced a 
maximum likelihood approach that calculates the cumulative probability of observing an ensemble 
of galaxies in a magnitude- limited survey for given a and M* . The probability of observing a galaxy 
of Mi(rrii,Zi) in a magnitude- limited redshift survey given a Schechter luminosity function is 

lQOA(M t -Mi(mi,Zi))(l+a) . ex pf_lQOA{M t -M i {m i ,Zi))\ 
Pi(Mi;M*) = @ ^ - x lo0 . 4(Mt _ M(m ^ ))(1+Q) . exp (_ 10 0.4(M,-M(m,*))W m ' ^ 

where 0j is the fraction of the angular area that is sensitive enough to detect the galaxy nii at 
redshift Zi and [m m i n , m max ] defines the magnitude sensitivity range of the survey. If the luminosity 
function remains the same for galaxies with different clustering properties at different redshifts, then 
the likelihood of obtaining an ensemble of N galaxies in the survey is 

i=N 

JC 4> =1[ Pi(M f ,M*). (10) 
i=i 

The faint-end slope a and characteristic magnitude M* in equation (8) are determined by maxi- 
mizing the likelihood function C^. 

The normalization fa of the Schechter luminosity function is determined separately, according 

to 

i=N I roc 

<t>* = Ys V ~ 1 / HM;M*)dM, (11) 

i= l / JM lim (zi) 

where V is the co-moving volume spcinncd. from ^min 

to 

z ma,x in the survey, <§>(M\ -M*) is the galaxy 
luminosity function defined in equation (8) with fa = 1, and Mn m (zi) is the corresponding absolute 
magnitude limit of the survey at the galaxy redshift Z{. The uncertainty in fa may be estimated 
by bootstrap resampling, that is, randomly re-sampling the same number of galaxies from the 
parent catalog a large number of times (allowing for duplications) and measuring the variation of 
fa between different re-sampled catalogs. 



3.3. The Stepwise Maximum Likelihood Approach 

A different maximum likelihood approach was developed by Efstathiou, Ellis, & Peterson 
(1988), one which does not require a functional form in the analysis and is similarly insensitive to 
the presence of large-scale density fluctuations in the survey. In this approach, the galaxy luminosity 
function is parameterized according to 



$(M) = fa(M k ), if M fc - ±AM <M<M k + ±AM, for k = 1, .., L. 



(12) 
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The probability of observing a galaxy of Mj(mj,£j), for which — ^AM < Mj(mj, Zj) < Mj, + 
±AM is, 

/fc=fc(m max ) 
^ <f) k (m, Zi) Am, (13) 
fc=fc(m min ) 

where fe(m m i n , Zj) and /c(m max , 2i) indicate the minimum and maximum absolute magnitude inter- 
vals defined by the apparent magnitude sensitivity range of the survey at Zi. The likelihood of 
observing N galaxies in a magnitude-limited survey is then 

i=N 

£^=11 Piimuzi). (14) 
i=i 

The best-fit 0fc's are obtained by maximizing the likelihood function defined in equation (14). 

In the simplest case where k(m m i n ) = ko and fc(m max ) = kf for i = 1-N, it is straightforward 
to show that oc m/%, where is the number of galaxies in each magnitude interval. Therefore, 
in the absence of large-scale density fluctuations, the SWML approach and the 1/V^ max approach 
are essentially the same. 



3.4. Systematic Uncertainties 

Photometric redshift techniques clearly provide an efficient means of identifying distant, faint 
galaxies in deep, wide-field surveys. It is, however, important to understand how the redshift 
measurement errors affect the derived luminosity functions of distant galaxies. Photometric redshift 
uncertainties propagate through the uncertainties in luminosity distances and consequently yield 
large uncertainties in the absolute magnitudes of individual galaxies. Here we neglect the effect 
of purely photometric uncertainties, because according to equation (4) a redshift uncertainty of 
5z = 0.1 will result in a magnitude uncertainty of AM ~ 0.4 at z ~ 1 due to the uncertainty in 
the luminosity distance. It is clear that photometric redshift errors dominate the uncertainties in 
the absolute magnitude calculations of individual galaxies. 

We performed a Monte Carlo simulation in order to understand the systematic uncertainties 
in the galaxy luminosity function calculations due to errors in photometric redshifts. First, we 
generated a model catalog of 8000 galaxies of different brightness at random redshifts between 
z = 0.1 and z = 1 according to an input Schechter luminosity function. Next, we assumed a redshift 
error function parameterized as a Gaussian distribution function of 1 a width 0.15 (1+Zi) and formed 
an "observed" redshift catalog by perturbing the input galaxy redshift Zi within the redshift error 
function. Finally, we determined the luminosity function for the galaxies at 0.5 < z < 0.8 using 
the 1/Vmax and STY methods. 

Figure 8 shows the derived galaxy luminosity function based on the simulated catalogs using 
different approaches. The input galaxy luminosity function is shown as the dash-dotted curve with 
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the selected M* and a indicated as a star in the inset. The solid circles represent measurements using 
the "observed" galaxy redshifts based on the 1/V max approach. The error bars are the associated 
1 a uncertainties estimated with a bootstrap re-sampling technique that takes into account both the 
sampling and redshift uncertainties. The dotted curve represents the best-fit Schechter luminosity 
function to the "observed" redshift catalog. The best-fit M* and a correspond to the upper solid 
point in the inset with the 99% uncertainties indicated by the dotted contour. 

It is clear from Figure 8 that large redshift errors together with the steep slope at the bright end 
of the galaxy luminosity function tend to flatten the observed luminosity function and result here in 
a best-fit M* that is 0.8 mag brighter than the input value. This confirms the statements above that 
more intrinsically fainter galaxies are scattered into the brighter end of the luminosity function than 
are intrinsically brighter galaxies scattered into the fainter end. This is equivalent to convolving 
the galaxy luminosity function with an error function in the galaxy absolute magnitudes induced 
by the redshift uncertainties. Evidently, redshift measurement uncertainties must be accounted for 
in order to obtain an accurate estimate of the intrinsic galaxy luminosity function. 



3.5. A Modified Maximum Likelihood Analysis 

In order to reduce the aforementioned bias we modified the likelihood analysis to incorporate 
the error functions of photometric redshift measurements into the probability calculations. Because 
objects are scattered in redshift space due to the relatively large redshift measurement errors, it is 
necessary to rewrite equations (9) and (13). The probability of observing a galaxy of Mj(mj, zi) in 
a magnitude-limited photometric redshift survey is now the cumulative probability for a galaxy of 
apparent magnitude rrii at z' to be identified at Zi according to the redshift error function pi(zi—z; Zi) 
defined in equation (3). Namely, it is a convolution of equations (9) and (13), respectively, with 
the error function of the photometric redshift measurement of the galaxy. Equation (9) defined for 
the STY approach is therefore 

J*f 10 0.4(Af.-M i (m i) z'))(l+a) . exp (_ 10 0.4 (M*-M,(m^'))) . - z > ;Zi ) dz > 



Pi (mi, zf, M* 



• f mmax r n 2/ l0OA(M t -M(m,z'))(l+a) . eX p(-10 - 4 {M.-M{m,z'))\ . p .U. _ z * . z .\ dz > dm ' 

(15) 

where Zf is the redshift limit of the grid search in the photometric redshift analysis, and equation 
(13) defined for the SWML approach is 

, , J Z/ <j> k (mi,z , )-p i (z i -z , ;z i )dz' 

^m^Zi)- . [io) 

T, k=k \ mini J V Mm, z>) ■ Pl ( Zi - z>; zi) dz'Am 

In cases where an accurate and precise spectroscopic redshift measurement is available, Pi(zi — 
z; Zi) = S(z — Zi) and equations (15) and (16) reduce to equations (9) and (13). The parameters 
that determine the Schechter luminosity function, a and M* in the STY approach and the 0fc's in 
the SWML approach, are then evaluated with the maximum likelihood functions (equations 10 and 
14). 
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We applied this modified likelihood analysis to the simulated catalog described in § 3.4 in order 
to determine how accurately we recover the intrinsic luminosity function with this prescription for 
the "empirical" error functions of the photometric redshifts. Here the redshift error function of 



results are shown in Figure 8 as well. The solid curve shows the best-fit Schechter luminosity 
function. The best-fit M* and a are the lower solid point in the inset with the 99% uncertainties 
indicated by the solid contour. The step function shows the best-fit stepwise luminosity function 
with the vertical bars indicating the associated 3 a, one-parameter uncertainties. Comparison of the 
simulation results for the luminosity function measurements with the input model shows that we 
are able to significantly reduce the systematic uncertainties (by ~ 0.7 mag) by taking into account 
an empirical redshift error function for each galaxy in the STY and SWML approaches. Although 
a lack of data means we could not constrain the faint-end galaxy luminosity function very well, 
we were able to recover the bright-end galaxy luminosity function to approximately one magnitude 
fainter than M*. 

3.6. The Galaxy Luminosity Functions from the LCIR Photometric Redshift Survey 

We applied the modified maximum likelihood analysis to determine the rest-frame i?-band 
galaxy luminosity function for ~ 3000 //-band selected galaxies identified at 0.5 < z < 1.5 in 
the LCIR survey. We adopted the redshift error function as formulated in equation (3), which 
is a convolution of the redshift likelihood function and a Gaussian kernel of width <r 2 (l + z) with 
o z = 0.08 as determined empirically from the comparison of photometric and spectroscopic redshifts 
described in § 2.3. 

The results of the analysis are presented in Figures 9 and 10 for the //-band detected galaxies. 
We calculated the luminosity functions for the total //-band selected sample (left panels) and for 
galaxies that have a best-fit spectral template of either E/S0 or Sab in the photometric redshift 
likelihood analysis (right panels). The solid curves in Figure 9 represent the best-fit Schechter 
luminosity functions at 0.5 < z < 0.75 (top panels), 0.75 < z < 1.0 (middle panels), and 1.0 < 
z < 1.5 (bottom panels) determined with the modified STY approach. Only galaxies identified in 
the HDFS region were included in the luminosity function calculations at 0.5 < z < 0.75 for the 
reasons discussed in § 2.4. For comparison, we include the r*-band luminosity function determined 
for galaxies at z < 0.2 from the Sloan Digital Sky Survey (SDSS; Blanton 2001; dotted curve) in the 
top-left panel. In addition, the best-fit luminosity functions of the LCIRs galaxies at 0.5 < z < 0.75 
are also indicated in the middle and bottom panels as the short-dashed curves. Because of the lack of 
sensitivity to intrinsically faint galaxies in the //-band survey, the solid curves shown in the middle 
and bottom panels were determined by fixing a to the best-fit value at 0.5 < z < 0.75 and allowing 
only Mft t to vary in the likelihood analysis. We have, however, attempted to determine both a 
and Mr. for galaxies at 0.75 < z < 1.0 and the results are shown as the long-dashed curves in the 
middle panels. The best-fit galaxy luminosity functions determined based on the modified SWML 



each galaxy is pi(z — zf, Zi) = 




f (1 + Zi) 2 } with a' z = 0.15. The 
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approach are shown as solid points in Figure 9. The vertical bars indicate the corresponding 1 a 
one-parameter uncertainties determined according to 2 In = 2 In £0 jmax — A% 2 with Ax 2 = 1. 

The good agreement at the bright end between the two best-fit Schechter luminosity functions 
(the solid and long-dashed curves) in each of the two middle panels of Figure 9 shows that despite 
the lack of sensitivity at the faint end, we can determine the bright-end galaxy luminosity function 
and therefore Mr^ fairly reliably using the iJ-band selected LCIR survey sample. In addition, 
a comparison of the solid and short-dashed curves in each panel indicates that there is little or 
no evolution in the rest-frame .R-band galaxy luminosity function for the entire ii"-band selected 
sample. The 99% error contours of the best-fit Mr„ and a using the modified STY analysis arc 
shown in Figure 10 for galaxies in different redshift intervals. It is clear from the top panels that the 
faint-end slope for galaxies at these redshifts is poorly constrained because our survey is insensitive 
to intrinsically faint galaxies at z > 0.75. 

We summarize the results of the modified STY analysis in Table 1, in which we list the 
redshift range, the best-fit faint-end slope a, Mr, — 5 log h, and 0* together with the associated 
1 a, one-parameter errors. The number of galaxies that are included in the analysis for all H- 
band selected galaxies and the if-band selected early-type galaxies are also listed. The results of 
the modified SWML analysis are summarized in Table 2, in which we list the rest-frame i?-band 
magnitude, Mr — 5 log h, and the best-fit log 4>(Mr) / (h 3 Mpc~ 3 ) together with the associated la, 
one-parameter errors in three different redshift ranges for the ii-band selected total and early-type 
samples. 

4. Discussion 

We have developed a technique to obtain a robust estimate of the intrinsic luminosity function 
using photometric redshifts. We adopt an "empirical" estimate of the redshift error function for each 
galaxy, which is a convolution of the redshift likelihood function that characterizes the uncertainty 
due to photometric errors and a Gaussian kernel of l-a width a z that characterizes the uncertainty 
due to template mismatch. This technique takes into account the non-gaussian error characteristics 
of photometric redshifts, and thus differs from the approach of Subbarao et al. (1996), who used 
a simple Gaussian to model the photometric redshift uncertainties. The results of the simulation 
described in § 3.4 show that our approach allows us to completely recover the galaxy luminosity 
function to the sensitivity limit of the survey. Here we compare the derived galaxy luminosity 
functions at different redshifts and discuss the implications for the evolution of galaxies since 
z ~ 1.5. 
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4.1. Galaxy Evolution at 0.5 < z < 1.5 

It is clear from Figure 9 that the depth of the ii-band survey is not sufficiently sensitive to con- 
strain the statistical properties of galaxies fainter than ~ at z > 0.75, and we are therefore unable 
to determine the evolution of the faint galaxy population based on the ii-band survey. However, 
the data are sufficient to determine the evolution of the bright galaxy population at 0.5 < z < 1.5 
based on the ii-band survey data. The evolution of this population has significant implications for 
understanding how galaxies form because different formation scenarios make distinct predictions for 
the number density evolution of massive galaxies. For example, under monolithic collapse scenarios 
(e.g. Eggen, Lynden-Bell, & Sandage 1962), massive galaxies form early on a dynamical timescale 
and have a constant co-moving space density and gradually declining intrinsic luminosities. In 
contrast, hierarchical formation scenarios (e.g. White & Rees 1978) predict that massive galaxies 
form through the merger of lower-mass galaxies over a Hubble time and therefore the co-moving 
space density of massive galaxies increases as the universe evolves. 

The luminosity functions presented in Figure 9 and Tables 1 and 2 show several interesting 
results. First, the best-fit Schechter luminosity functions agree with the stepwise luminosity func- 
tions to within the measurement uncertainties at all redshifts, including the luminosity function at 
0.5 < z < 0.75 where our survey is sufficiently sensitive to identify galaxies as faint as ~ 0.2 L*. 
Second, there exists little or no evolution either in Mr^ or in for the total ii-band selected 
population at 0.5 < z < 1.5. Third, we find a moderate luminosity evolution for galaxies with 
an SED best characterized as E/S0 or Sab galaxies. Specifically, our measurements obtained from 
the modified STY approach suggest that from z ~ 1 to z ~ 0.5 these early-type galaxies have 
faded by ~ 0.5 mag. This is consistent with the expected luminosity evolution for a galaxy formed 
at Zf = 30 and following an exponentially declining star formation rate model with r = 1 Gyr 
(McCarthy et al. 2001; Chen et al. 2002). Finally, we find only a mild evolution in 0* between 
z w 0.5 and z ~ 1 (~ 40% decrease with redshift) for the color-selected early-type galaxies, but we 
cannot argue against the no evolution scenario at more than the 2 a level. When we consider the 
entire ii-band selected sample at 0.5 < z < 1.5, our estimates of (ft* are completely consistent with 
a no evolution scenario. 

It appears that neither the total ii-band selected sample nor the color-selected early-type 
galaxies have evolved significantly in intrinsic luminosity or space density since z ~ 1.5, contrary 
to the predictions of the hierarchical formation scenarios. But as with all luminosity function 
measurements, it is clear that there exists a strong degeneracy between Mr„ and 0* at z > 0.75 in 
our analysis, because the ii-band survey is not sensitive to sub-L* galaxies at this redshift range 
and our estimates for these two parameters were determined from the space density of intrinsically 
bright galaxies with a fixed faint-end slope a. For example, a bright Mji t together with a small 0* 
can produce the same bright-end luminosity function as a faint Mr^ with a large </>*. To estimate 
how significantly early-type galaxies have evolved since z ~ 1.5, we therefore calculated the rest- 
frame co-moving i?-band luminosity density £r for galaxies brighter than M m - m using the best-fit 
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Schechter luminosity function, 

/oo 
L-<f>(L;L*)d(L/L*) (17) 
—"tnin 

(0.4 In 10) ■ 0* ■ L* • W V-HM*-M)(2+a) . exp (_ 10 0.4(M,-M) )M (lg) 

-oo 

This represents an integrated quantity that characterizes the brightest galaxy population and is 
least sensitive to faint galaxies. 

We test different formation scenarios using the estimated co-moving luminosity density evolu- 
tion of the bright, red galaxies. Including the luminosity function estimated for early-type galaxies 
at z = 0.3 in the CNOC2 survey (Lin et al. 1999), we first compared £r calculated for a constant 
M m i n at all redshifts. The results are presented in the left panel of Figure 11, which shows a flat £r 
versus z for M min = —20.5 (squares) , —21.0 (triangles), and —21.5 (circles). These roughly corre- 
sponds to L*, 1.6 L*, and 2.5 L*, respectively, at z ~ 0.3. However, because stars evolve with time, 
the corresponding luminosity evolution must be accounted for. To distinguish between a passively 
evolving population and hierarchical mergers, we estimated the amount of brightening (A Mr) with 
increasing redshift (or with decreasing age of the universe) for different stellar evolution models 
(see Table 3 for a summary). We considered a wide range of stellar evolution scenarios, from a 
single burst at Zf = 30, to an exponentially declining star formation rate model with Zf = 5 and 
r = 1 Gyr. The former represents a classic passive evolution mode under the monolithic collapse 
scenariol. The latter represents the closest resemblance of a more continuous star-forming scenario 
in which some residual star formation is taking place at z < 3, while there is still enough time for 
the red galaxy population to be present by z ~ 1. We then calculated £r for a suite of evolving 
M m i n (z) according to different stellar evolution recipes. 

The goal of this analysis is to examine whether or not mergers are the dominant process for the 
formation of red galaxies. Mergers would lead to brightening with time (smaller systems combined 
to form bigger ones), which compensates the expected fading due to stellar evolution. We would 
therefore expect to observe an increasing £r with decreasing redshift using (A Mr) presented in 
Table 3, if these red galaxies formed through a sequence of mergers. A flat redshift distribution of 
£r would then imply that mergers are insignificant in the redshift interval tested. 

We present the results in Figure 11. The middle panel is for the Zf = 30 single burst scenario 
and the right panel is for the Zf = 5 exponentially declining star formation rate model. The results 
in these two panels have been scaled to have consistent M m [ n at z = 0.3 as in the left panel. 
In addition, we exclude the brightest subsample (circles in the left panel) from these two panels, 
because the number of galaxies at the very bright end of the luminosity functions is very small 
beyond z ~ 0.75. Specifically, there are only eight galaxies in the red sample that are brighter than 
M R = -23 at 1.0 < z < 1.5. 

Figure 11 shows that £r evolves relatively slowly with redshift for different bright subsamples, 
and that it could have increased by at most a factor of « 6 from z ~ 1.2 to z ~ 0.3 for the 
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brightest early-type sample under the most extreme stellar evolution model presented in the right 
panel. Taking into account the fact that these models do not have dust included, which may further 
reduce the amount of fading in intrinsic luminosity with redshift, we find no evidence to support 
that mergers dominate the formation of early-type galaxies over 0.5 < z < 1.5. This is consistent 
with the findings of Firth et al. (2002), who directly compared the observed abundances of red 
galaxies at z ~ 1 with those from different semi-analytical models. They found that the models 
underpredict the number of red galaxies observed in the LCIR survey by at least a factor of three. 

The rest-frame co-moving i?-band luminosity densities for all if -band selected galaxies at 
different redshifts are presented in Figure 12. As in Figure 11, we calculated £r for three constant 
thresholds, M m i n = —20.5 (squares) , —21.0 (triangles), and —21.5 (circles). Because galaxies in 
the total -ff-band selected sample have a wide range of star formation history, we cannot remove 
the effects of pure luminosity evolution by adopting a simple M m { n (z). The results of Figure 12 
show that there exists only a mild evolution in the rest-frame co-moving i?-band luminosity density 
for the entire galaxy population over 0.5 < z < 1.5. 

It is also important to note that our analysis is based on an i^-band selected galaxy sample, 
which roughly corresponds to galaxies selected in the rest-frame J, /, and R bands at z ~ 0.5, 
0.9, and 1.3, respectively. A potential systematic bias in our measurements is clearly the bandpass 
selection effect which, when combined with variation in star formation rate and stellar evolution 
with time, makes it difficult to interpret our results, particularly for the total ff-band selected sam- 
ple. At z ~ 1.3, the H-hand survey preferentially selects optically bright (and therefore younger) 
galaxies, while it tends to select near-infrared bright (and therefore older) galaxies at z ~ 0.5, when 
the universe was twice as old. Because the total H-band selected sample consists of galaxies of 
different star formation histories in different redshift intervals, it is possible that galaxies of differ- 
ent colors, selected at different rest-frame wavelengths at different redshifts conspire to produce a 
luminosity function that appears to show little evolution. Consequently, it is impossible to rule out 
any evolution scenarios without a detailed understanding of the intrinsic properties of the galaxies 
in the total i^-band sample. On the other hand, the bandpass selection bias has a much less serious 
effect for the early-type galaxy sample. Only galaxies with an observed SED (spanning a spectral 
range from the V band or U for the HDFS objects through the H band) that is consistent or 
redder than a typical Sab galaxy were included in the analysis at all redshifts. Our measurements 
represent a true, empirical measure of how "red" galaxies evolve with time, and the results of our 
analysis suggest that early- type galaxies have not evolved significantly since z ~ 1.5. 



4.2. Comparison with Previous Studies 

The luminosity functions presented in § 3.6 cover the redshift range 0.5 < z < 1.5. The upper 
end of this range extends beyond the reach of traditional spectroscopic redshift surveys, the most 
ambitious of which stretch to z = 1.2 but become significantly incomplete beyond z = 1 (Lilly et al. 
1995; Cowie et al. 1996; Cohen 2002). Various groups have measured evolution in the luminosity 
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function using spectroscopic redshifts over the redshift range < z < 1 for galaxies of different 
colors (e.g. Lilly et al. 1995), different spectral types (e.g. Heyl et al. 1997; Cohen 2002), or different 
morphology (Im et al. 1999). In addition, some groups have attempted to measure the luminosity 
function using photometric redshifts for galaxies at z < 4 (Gwyn & Hartwick 1996; Mobasher 
et al. 1996; Connolly et al. 1997; Sawicki et al. 1997). Different sample selection criteria make 
a quantitative comparison between these measurements difficult. Here we present a qualitative 
comparison between our measurements and the previous ones. 

First, we consider previous measurements that were conducted in the rest-frame R band. The 
local luminosity density is now well constrained by both the 2dF Redshift Survey (Cross et al. 2001) 
and the SDSS (Blanton 2001). It is encouraging that these surveys yielded consistent results in the 
B band. The i?-band luminosity density available from the SDSS is most directly relevant to our 
results from the -ff-selected LCIR survey, and is included in Figure 12 for the luminosity density 
evolution of the total .ff-band selected sample. If we consider only the LCIR measurements at 
0.5 < z < 1 and the SDSS measurement at z < 0.2, then the evolution in the luminosity density 
from z ~ to z ~ 1 is less than 0.2 dex for any reasonable choice of the faint-end slope, a, indicating 
a very slow evolution. 

Next, we compare our results with those from a number of very deep redshift surveys that 
covered a relatively smaller area but were sensitive to galaxies at redshift beyond z ~ 1. These are 
the Canada- France Redshift Survey (CFRS; e.g. Lilly et al. 1996), a deep ET-band survey by Cowie, 
Songaila, & Barger (1999), and the Caltech Faint Galaxy Redshift Survey (CFGRS; Cohen et al. 
2002), all of which showed that the bulk of the evolution in the luminosity function occured in bluer 
spectral types. In addition, they showed that galaxies with redder SEDs evolved more slowly, if at 
all, over this redshift range, a result with which the LCIR survey is at least qualitatively consistent. 
Note however that these results may change if adopting a different color selection criterion for red 
galaxies at high redshifts (Kauffmann, Chariot, & White 1996) or the incompleteness correction 
for absorption feature dominated galaxies was incorrect (Totani & Yoshii 1997), two effects that 
work in opposite ways. 

A more quantitative comparison is possible with Lilly et al. (1996), which gave the evolution of 
the luminosity density between z = and z = 1 for different CFRS samples. In particular, Figure 
12 is most directly comparable to Figure 1 in Lilly et al. (1996). The CFRS luminosity densities 
were computed at rest-frame 2800 A, 4400 A, and 1 fim by interpolating their BVIK photometry, 
which showed a decline by approximately 0.5 dex at 4400 A and 0.4 dex at 1 fim from z ~ 1 to 
z ~ 0.* This gives A log t /A log(l + z) = 1.7 ± 0.5 at 4400 A and 1.3 ± 0.4 at 1 fim. We find from 
the LCIR survey that Alog^ /Alog(l + z) = 0.6 ± 0.1 at 6800 A, if we assume a faint-end slope 
a = —1.0. If we assume a very steep faint-end slope a = —1.5, then we measure still only 1.0 ±0.2. 
Although our results are marginally consistent with the CFRS results, the i^-selected LCIR survey 



*We have subtracted 0.2 dex from the original measurements in Lilly et al. to account for the differences in the 
adopted cosmological models. 



- 22 - 



appears to favor slower evolution in the overall galaxy population at z < 1 (at the « 1 a level). 

The slow evolution measured here is similar to the /-band evolution measured by Cowie, 
Songaila, & Barger (1999) from K-selected spectroscopic surveys at z < 1.2. Based on a much 
larger sample covering a broader redshift range, we find that this slower evolution at red wavelengths 
extends to at least z = 1.5. Our results are also consistent with the interpretation of the luminosity 
function evolution derived from the CFGRS, which represents by far the largest and most complete 
spectroscopic sample of galaxies at z > 0.8. It is, however, still limited by small sample size (144 
galaxies at 0.8 < z < 1.05 and 18 galaxies at 1.05 < z < 1.5). This illustrates the difficulties of 
carrying out a complete spectroscopic survey over a large area to very faint magnitude limits. 

Next, we compare our results with those estimated based on photometric redshifts. At 
1 < z < 2, photometric redshift surveys prevail but have so far been limited to small areas of 
the sky. In addition, it is important to include near-infrared photometry in order for accurate pho- 
tometric redshift measurements to be plausible. Connolly et al. (1997) showed that uncertainties 
of photometric redshifts are reduced by 40%, when incorporating near-infrared photometry in the 
photometric redshift analysis. Indeed, the reddest galaxies known in this reshift interval (Elston 
et al. 1988; McCarthy et al. 1992; Hu & Ridgeway 1994) are often exceedingly faint in all optical 
bandpasses bluer than the R band, making photometric redshift measurements very uncertain. The 
luminosity function from Connolly et al. was for a J-band selected galaxy sample over 1 < z < 2 
in the Hubble Deep Field North (HDFN). Their luminosity function was computed, however, at 
rest-frame 2800 A, which is more sensitive to young stellar populations than the rest-frame R band. 
Furthermore, over the redshift range probed by the LCIR survey (0.5 < z < 1.5), the HDFN is pri- 
marily sensitive to the faint end of the luminosity function, while the LCIR survey is sensitive only 
to the bright end. Finally, the total survey area of HDFN was only 4.48 arcmin 2 , corresponding to a 
comoving volume of ~ 1.3 x 10 4 Mpc 3 over the 1 < z < 2 interval. In contrast, our analysis is based 
on a survey volume that is more than two orders of magnitude larger than in the HDFN. Therefore, 
we have sufficient numbers of bright galaxies to z = 1.5 to probe their evolution in greater detail. 
The larger volume is necessary to accurately measure the mean space density because the bright, 
red galaxies are now known to be strongly clustered (Daddi et al. 2000; McCarthy et al. 2001; Firth 
et al. 2002). 

Other groups have attempted to study the evolution of the galaxy luminosity function at z < 1 
based on a deep photometric redshift survey that incorporates near-infrared photometry (Fried et 
al. 2001; Poli et al. 2001), but these galaxy samples were selected in the / band, which lends a 
bluer tint to the overall mix of SEDs than the .ff-selected LCIR survey galaxies. It is therefore not 
straightforward to compare our results with these other photometric redshift surveys. 

Finally, we consider measurements obtained for early-type galaxies selected on the basis of 
quantitative morphological criteria such as smoothness, symmetry, and bulge-to-total ratio. Im et 
al. (1999) selected a sample of 145 E/S0 galaxies from deep HST imaging in the Groth Strip (Groth 
et al. 1994) Using a subset with spectroscopic redshifts as calibrators, they estimated redshifts for 
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the remaining galaxies using the V — I color alone. The resulting luminosity functions (measured 
in the rest-frame B band) suggest slow evolution in the E/SO population at z < 1, in qualitative 
agreement with the results discussed earlier. Brinchmann et al. (1998) and Menanteau et al. 
(1999) have, however, suggested that the bulk of morphologically selected early-type galaxies at 
z < 1 have bluer colors than those predicted by different stellar evolution models, but the results 
are inconclusive because they depend sensitively on the adopted models and how photometric 
uncertainties are incorporated in these models. 

In summary, we find that the rest-frame, co-moving optical luminosity density of the entire 
galaxy population has evolved only moderately (declining by a factor of « 2) from z ~ 1 to z ~ 0. 
The rest-frame, co-moving optical luminosity density of the color-selected, early-type galaxies has 
not evolved significantly since z ~ 1. Our new measurements of luminosity density evolution is 
based on a significantly larger sample selected in the near-infrared than previous surveys selected 
at a range of wavelengths. It is consistent with previous findings, but we place a much tighter 
constraint on the apparant lack of evolution in the co-moving luminosity density of early-type 
galaxies. 

5. Summary 

We have obtained an estimate of the galaxy luminosity function at z > 1 using photometric 
redshifts for galaxies identified in the HDFS and CDFS regions of the LCIR survey. We have also 
calculated the luminosity functions at redshifts between z ~ 0.5 and z ~ 1.0 using the same galaxy 
sample selected uniformly at near- infrared wavelengths. The luminosity function analysis is based 
on ~ 3000 i^-band selected galaxies with apparent magnitude 17 < H < 20 at 0.5 < z < 1.5. 
We demonstrate that our photometric redshift measurements are accurate with an RMS dispersion 
between the photometric and spectroscopic redshifts of a z j (1+z) ~ 0.08. In addition, we show that 
the systematic uncertainty inherent in the luminosity function measurements due to uncertainties in 
photometric redshifts is non-negligible and therefore must be accounted for. We develop a technique 
to incorporate photometric redshift error functions of individual galaxies in the estimates of the 
intrinsic luminosity function. We show that this technique allows us to completely recover the galaxy 
luminosity function to the sensitivity limit of the survey and therefore to study the luminosity and 
space density evolution of near-infrared selected galaxies over a large redshift interval. The results 
of our study are: 

1. On the basis of galaxies identified in the Hubble Deep Field South and Chandra Deep Field 
South regions, we find that (7.3 ± 0.2) % and (16.7 ± 0.4) %, respectively, of the iJ-band detected 
galaxies of H < 20 are at z > 1. 

2. The galaxy luminosity functions determined via a modified SWML method are consistent 
with the Schechter functions determined via the modified STY method. 

3. The best-fit Mr, and 0* for the total .ff-band selected population are consistent with a 
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no-evolution model at 0.5 < z < 1.5. 

4. The evolution of the co-moving luminosity density of the .ff-band selected galaxies is char- 
acterized by Alog£/Alog(l + z) = 0.7 ± 0.2 at rest-frame 6800 A. This is marginally consistent 
with the CFRS results, but favors a slower evolution (at the « la level). 

5. Galaxies with spectral energy distributions best characterized as E/S0 or Sab appear to 
have faded by ~ 0.5 mag from z ~ 1 to z ~ 0.5, consistent with the expected luminosity evolution 
from an exponentially declining star formation rate model. In addition, the space density of the 
these color-selected early-type galaxies has evolved only mildly 1 (~ 40% decrease with 
redshift). 

6. The derived rest-frame co- moving i?-band luminosity density £r of color-selected early- type 
galaxies exhibits only moderate evolution with redshift. Specifically, we find that under the most 
extreme stellar evolution scenario the £r of early- type galaxies brighter than L* (1.6 L*) could have 
increased by at most a factor of « 3(6) from z ~ 1.2 to z ~ 0.3. Our results place an upper limit 
to the merger rate of bright galaxies over this time interval. 
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Table 1. Results of the modified STY Analysis 1 



All E/SO + Sab 



Mr, fa Number Mr, fa Number 

Redshift Range a —5 log h (/i 3 Mpc -3 ) of Galaxies a —5 log h (/i 3 Mpc -3 ) of Galaxies 



[0.50,0.75] -l.OOlgoa -21.03^ 0.0142±g-ggi§ 978 -0.20j#o4 -20.431°^ 0.0127±g-ggi| 606 

[0.75,1.00] -0.50l° : °« -21.03±g : °t 0.0142±g : ggi| 1239 +0.20^™ -20.83$$ 0.0070±g : g°g| 752 

-i.ooi8:88 -2i.44t8:8? o.oiosi^ 1239 -0.2018:88 -2i.o4±g:gi 0.00721s 752 

[1.00,1.50] -1.0018:88 -20.9218:8^ 0.027118:88^ 791 -0.2018:88 -20.6718^ o.omi™^ 4 58 



a Numbers with zero errors indicate that the parameters were fixed at the listed values in the likelihood analysis. 
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Table 2. Modified SWML Galaxy Luminosity Function, log 4>{M R ) / '(/i 3 Mpc~ 3 ) 







All 






E/S0 + Sab 




Mr — 5 log h 


[0.50,0.75] 


[0.75,1.00] 


[1.00,1.50] 


[0.50,0.75] 


[0.75, 1.00] 


[1.00,1.50] 


-23.23 
-22.73 
-22.23 
-21.73 
-21.23 
-20.73 
-20.23 
-19.73 


-5 52 +L07 
°- oz -0.01 

'yoH _ 0.32 
O.lO_ 010 

-2 Q6 +0 - 18 
z - yo -0.06 

9 97+O.O9 

z - z ' -0.03 

1 OQ+0.02 
i - yo -0.06 

_1 Q4+0.02 
i - yq: -0.02 
9 91 +0.03 
z - zi -0.06 


_ 7 ic+1.42 
' • lo -0.16 

-5 m+ - 65 

°- ui -0.58 
90+0.28 
°- zu -0.02 
9 97+O.OI 
z - z ' -0.09 
9 9-1+0.04 
z - zl -0.02 
9 1 c:+0.03 
z - io -0.05 


-4 fi3+ - 77 

4 - D,3 -0.04 

_2.99i a02 
z - 04 -o.oi 


c 90 + I.2O 
J - zo -0.01 
_o 7 n+0.48 
°-' u -0.04 

Z - DO -0.02 

-2 1Q+ 001 
z - iy -0.16 

9 -, o+0. 03 
- z -18_0.01 
_n tro+0.03 

z - oo -0.07 


_4 ^7+0.22 

4 7^+0-74 
4 -'°-0.97 
9 oq+0.03 
- Z -88_0.14 

-2 64 +a04 
z - D4 -0.05 

9 7-1+0.04 

z -' i -0.04 

9 9 c:+0.06 

z - zo -0.04 


_ 4 O9+0.83 
4 - yz -0.15 

, 4fi +o.oi 
°-^ D -o.oi 

-2 60 +0 - 04 

z - ou -0.03 



Table 3. Predicted Brightening (A Mr) with Redshift for Different Star Formation Histories 



Scenario 


z f 


z = 0.3 


z = 0.625 


z = 0.875 


z = 1.25 


single burst with A t = 1 Gyr . . 


. 30 


-0.30 


-0.53 


-0.69 


-0.94 




10 


-0.32 


-0.55 


-0.74 


-0.98 


exp(— t/r) with r = 1 Gyr 


30 


-0.30 


-0.60 


-0.81 


-1.15 




10 


-0.34 


-0.64 


-0.88 


-1.24 




5 


-0.35 


-0.70 


-0.96 


-1.40 
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Fig. 1. — Median residuals of "optimal" (m opt ) vs. 4"-aperture photometry (m ap ) in the V, R, I, 
and H bandpasses for the if-band selected objects in the CDFS region. Horizontal bars show the 
magnitude bin size. Vertical bars mark the 1 a scatters of the residuals for individual magnitude 
bins. Straight line segment shows m opt = m ap . 
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Fig. 2. — Sample results of the redshift likelihood analysis for four spectroscopically identified 
galaxies in the HDFS region. The top panel of each pair shows the photometric measurements 
and uncertainties as solid points with error bars. The solid curves represent the best-fit spectral 
templates with the spectral type indicated in the upper-left corner. The open squares indicate 
the predicted fluxes in individual bandpasses, estimated by integrating the best-fit template over 
the corresponding system transmission functions. The dashed lines indicate the zero flux level. 
The bottom panel of each pair shows the redshift likelihood functions. The estimated photometric 
redshift (z p hot) and spectroscopic redshift (z S pcc) are indicated in the lower-right corner of each of 
the top panels. 
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Fig. 3. — Left: Redshift histogram of the .ff-band detected galaxies in the HDFS region. The 
hatched region indicates the redshift distribution of the I — H > 3 galaxies. Solid curves indicate 
the predicted redshift distribution based on the model that best fits the number-magnitude relation 
presented in Chen et al. (2002). Right: Cumulative redshift distributions of the fraction of the H- 
band selected galaxies with H < 20. The open and hatched histograms indicate galaxies identified 
in the HDFS and CDFS regions, respectively. We find that between ps 7% and 17% of the //-band 
selected galaxies are at z > 1. 



-34- 




Fig. 4. — Left: Comparisons of photometric (z p hot) an d spectroscopic (z spcc ) redshifts for the 
HDFS (solid points) and CDFS (open circles) objects. Right: The corresponding redshift residual 
distributions for the HDFS (open histogram) and CDFS (shaded histogram) objects. The solid 
line in the left panel indicates z p hot = %>ec- There are 67 objects in the HDFS region with known 
spectroscopic redshifts. We find that photometric redshifts are accurate to within a z /(l + z) w 0.08 
at z < 1. Based on 34 CDFS objects with reliable photometric and spectroscopic redshifts, we 
find that photometric redshifts estimated without U and B photometry are accurate to within 
a z /(l + z) « 0.08 at z spcc > 0.75. 
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Fig. 5. — Sample results of the Monte Carlo simulations to determine photometric redshift un- 
certainties for four different spectral types due to uncertainties in the galaxy photometry. Based 
on 1000 realizations of a galaxy of H = 20 at z = 1.2 in the HDFS region, the left panels show 
histograms of the residuals between the input galaxy redshift z- ul and the best-fit photometric red- 
shift 

^out; and the right panels show cumulative fractions of the simulated photometric redshifts 
versus absolute redshift residuals. In each realization, we perturbed the photometric measurements 
within the 1 a photometric uncertainties known for the individual optical and H images in the 
HDFS region. We find that for a galaxy of H = 20 at z = 1.2 in the HDFS region there is an 
f=s 10% chance that it would be misidentified as a starburst galaxy at z ~ 4, if the spectral shape 
mimicked an E/S0 galaxy. 
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Fig. 6. — Summaries of the Monte Carlo simulations to determine photometric redshift uncertainties 
due to uncertainties of galaxy photometry versus redshift and iJ-band magnitude. For a given input 
redshift and i7-band magnitude, each panel now shows the sum of the residual histograms over four 
spectral types as shown separately in Figure 5. The vertical axis is on a logarithmic scale as shown 
in Figure 5. The upper panels are results based on the filter combination and image sensitivities 
of the HDFS galaxies and the lower panels for the CDFS galaxies. 
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Fig. 7. — Cumulative fractions of the simulated photometric redshifts versus absolute redshift 
residuals for the four z = 1.2 cases presented in Figure 5. The left panels show the simulation 
results for the HDFS galaxies and the right panels for the CDFS galaxies. At z = 1.2, it is clear 
that photometric redshifts are accurate to within Az = 0.1 for the CDFS galaxies of H < 21, and 
that they become very uncertain for the HDFS galaxies at H ~ 21 because of their relatively larger 
photometric uncertainties. Specifically, there is « 30% chance that the photometric redshifts of 
these galaxies are off by more than Az = 0.1. 
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Fig. 8. — Simulation results showing systematic uncertainties in the luminosity function calcula- 
tions due to the relatively large redshift uncertainties. The dash-dotted curve indicates the input 
(intrinsic) galaxy luminosity function with the selected Schechter-function parameters, M* and a, 
indicated by the star in the inset. The solid circles represent measurements using the "observed" 
redshift catalog based on the l/V mSuX approach. Error bars are the associated la uncertainties 
estimated using a bootstrap re-sampling technique that takes into account both the sampling and 
redshift uncertainties. The dotted curves represent the best-fit Schechter luminosity function based 
on the STY approach without including the redshift error functions of individual galaxies in the 
likelihood analysis. The best-fit M* and a are indicated by the upper solid point in the inset 
with the dotted contour showing the 99% uncertainties. The best-fit Schechter luminosity function 
based on our modified maximum likelihood analysis is shown as the solid curve. The best-fit M* 
and a are indicated by the lower solid point in the inset with the solid contour showing the 99% 
uncertainties. The step function shows the results of the modified SWML method with the vertical 
bars indicating the associated 1 a, one-parameter uncertainties. 
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Fig. 9. — The rest-frame i?-band galaxy luminosity functions calculated using the modified STY 
(curves) and SWML (solid points) approaches. The measurements are for the total If -band selected 
galaxy sample (left panels) and for galaxies that have a best-fit spectral template of either an E/SO 
or an Sab galaxy in the photometric redshift likelihood analysis (right panels). The results shown 
in the top panels were obtained based on galaxies detected in the HDFS regions only, excluding the 
CDFS galaxies at z < 0.75 for which the photometric redshifts become unreliable due to the lack of 
U and B photometry. The horizontal bars indicate the bin size in absolute magnitude AMr. The 
vertical bars associated with the solid points represent the 1 a, one-parameter uncertainties in the 
SWML analysis. The solid and long-dashed curves represent, respectively, the best-fit Schechter 
luminosity functions with and without a fixed at the fiducial value at 0.5 < z < 0.75. The fiducial 
luminosity function at 0.5 < z < 0.75 for each sub-sample is shown as the short-dashed curve in the 
middle and bottom panels for comparison. For comparison, we also include in the top-left panel the 
r* luminosity function derived for galaxies at z < 0.2 from the SDSS (Blanton et al. 2001; dotted 
curve) . 
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Fig. 10. — The 99% error contours of the best-fit Mr, and a in the modified STY approach for 
the entire ff-band selected galaxy sample (left panels) and the //-band selected early-type galaxies 
(right panels) at 0.5 < z < 0.75 (bottom panels) and 0.75 < z < 1.0 (top panels). The cross in 
each panel indicates the best-fit Mr, and a for the total //-band sample at 0.5 < z < 0.75. 
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Fig. 11. — Redshift evolution of the rest-frame co- moving R-h&nd luminosity density Ir for color- 
selected early- type galaxies brighter than M m \ n . We show in the left panel the results calculated 
by adopting a constant M m j n versus redshift for three different threshold values: M m - m = —20.5 in 
squares, —21.0 in triangles, and —21.5 in circles. The z = 0.3 values were derived based on the 
galaxy luminosity function determined by Lin et al. (1999) from the CNOC2 early-type sample. We 
also calculated £r by adopting an evolving M m j n in order to remove the effect of pure luminosity 
evolution of stars. The middle panel shows the results for M m i n estimated under a single burst 
scenario for a galaxy formed at ~ 30. The right panel shows the results for M m - m estimated 
under a 1-Gyr exponentially declining star formation rate scenario for a galaxy formed at Zf ~ 5. 
The predicted brightening with redshift under different scenarios is listed in Table 3. The results in 
the two right panels have been scaled to have consistent M m i n at z = 0.3 as those in the left panel. 
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Fig. 12. — Redshift evolution of the rest-frame co- moving i?-band luminosity density for all .ff-band 
selected galaxies brighter than M m \ n . The results were calculated by adopting a constant M m - m 
versus redshift for three different threshold values: M m - m = —20.5 in squares, —21.0 in triangles, 
and —21.5 in circles. The z = values were derived based on the galaxy luminosity function 
determined by Blanton et al. (2001) from the SDSS sample. It shows that the luminosity density 
decreases from z ~ 1 to z ~ by approximately Alog^ /Alog(l + z) = 1.0 ± 0.2 at 6800 A. To 
examine how sensitive this evolutionary slope is to galaxies at the faint end, we extrapolated the 
calculation to M m j n = —17.5 by adopting the best-fit a = — 1.2 of the SDSS luminosity function at 
z = and a = — 1 (—1.5) for the LCIR luminosity functions at higher redshifts. The results are 
shown in open (dotted) circles and we find AlogZ /Alog(l + z) « 0.6-1.0. 



